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The  Dynamics  of  Formation  and  Evaporation  of 
Mixed  Alkali  Halide  Nanocrystals: 

A  Case  of  Comparable  Lattice  Energies 


Temer  S.  Ahmadi  and  Mostafa  A.  El-Sayed 
Department  of  Chemistry  and  Biochemistry 
University  of  California,  Los  Angeles,  California  90024 


Abstract 


The  relative  intensity  distribution  of  the  [RbnKi4-nIl3]”*‘  mixed 
nanocrystals  containing  a  "magic"  number  of  14  metal  cations  and  13 
iodide  anions  is  examined.  These  nanocrystals  were  generated  through 
sputtering  of  mixed  Rbl  /  KI  solids  using  fast  atom  bombardment,  and 
analyzed  by  use  of  a  double-focusing  sector-field  mass  spectrometer.  The 
relative  intensities  of  the  observed  peaks  is  found  to  have  a  'bimodal' 
distribution  with  a  maximum  near  ns?.  This  observation  is  explained  by 
assuming:  (a)  a  rate  of  formation  that  is  proportional  to  the  number  of 

O'les 


isomers  that  can  be  formed  for  each  mixed  nanocrystal,  and  (b)  a  rate  of 
dissociation  that  is  insensitive  to  the  composition  of  the  mixed  nanocrystal. 
The  later  assumption  leads  to  the  conclusion  that  the  internal  energies  of 
these  mixed  nanocrystals  are  similar  and  independent  of  their  composition. 
This  conclusion  is  supported  by  the  results  of  Monte  Carlo  simulation 
annealing  as  well  as  by  the  similarity  of  the  known  lattice  energies  of  KI 
and  Rbl  bulk  crystals. 


I.  Introduction 

Because  of  their  elementary  electronic  properties,  the  alkali  halide 
clusters  have  been  extensively  studied  both  theoretically  and 
experimentally In  1981,  Campana^*^  et.  al.  used  Xenon-ion 
bombardment  of  Csl  to  generate  cluster  ions  of  the  type  [MnXn-ll"''  (for  n 
values  up  to  70)  which  were  detected  by  a  secondary-ion  mass 
spectrometer.  They  noted  an  irregularity  in  the  relative  abundance  of  the 
cluster-ions,  indicating  small  deviations  from  the  expected^^’^^  monotonic 
decrease  in  ion  intensity  as  the  cluster  size  increased  (since  larger  clusters 
require  higher  number  of  collisions  for  their  nucleation).  These  deviations 
or  "local  maxima"  occur  at  n=:  14,  23,  38,  and  63  .  These  n  values  are  called 
the  "magic  numbers"  and  are  attributed  to  the  formation  of  relatively 
stable27'3 13X3X3,  3X3X5,  3X5X5,  and  5X5X5  cuboid  structures.  This 
relative  stability  diminishes  the  rate  of  their  decomposition  prior  to  their 
detection,  thus,  allowing  them  to  have  relatively  stronger  peak  intensities 
in  the  mass  spectrum. 
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A  direct  determination  of  the  structure  of  these  cluster-ions  has  not 
been  possible;  however,  Martin^ ^  used  total  energy  calculations  to 
determine  the  most  stable  structures  of  [NanCln-ll"^  cluster-ions  (n=  1-18), 
For  n  =14,  one  of  the  "magic  numbers",  he  found  a  3X3X3  face-centered- 
cubic  lattice  structure.  Calculations  by  others  have  confirmed  Martin's 
work^’^^.  The  mass  spectra  of  all  of  the  alkali  halides  have  the  same 
"magic  numbers",  indicating  a  face-centered-cubic  (fee)  structure  for  all  of 
them  irrespective  of  their  bulk  crystalline  structures,  e.g.  bulk  NaCl,  Nal, 
NaF,  KF,  KI,  and  KCl  crystallize  into  fee  structures,  whereas  bulk  CsCl  and 
Csl  form  a  body-centered-cubic  (bcc)  lattice^®’^^. 

The  structure  and  relative  yield  of  small  mixed  alkali  halide  clusters 
have  also  been  examined  Diefenbach  and  Martin found  that 

cluster-ions  of  mixed  alkali  halides  are  stable  even  for  compositions  that 
do  not  form  solid  solutions  in  the  condensed  phase  (e.g.  Nal  and  Csl). 
Using  the  polarizable-ion  model  ,  they  calculated  the  total  energies  and  the 
corresponding  geometries  of  these  clusters.  It  was  found  that  mixed 
clusters  with  ions  of  comparable  sizes(e.g.  Rb+,  Cs+)  give  the  most  stable 
structures'^. 

In  an  effort  to  examine  the  effect  of  defects  on  the  formation  and 

stability  of  nanocrystals,  we  have  previously^ ^  studied  the  addition  of  an 

Agl  molecule  to  the  Ag  cluster-ion.  The  structure  of  the  mixed  clusters 

n 

were  deduced  from  the  observed  evaporation  dynamics  and  the 
interactions  resulting  from  the  adhesion  of  the  Agl  molecules  to  the 
metallic  clusters.  In  the  present  work,  we  experimentally  determined  the 

intensity  distribution  of  mixed  alkali  halide  nanocrystals  as  Rbl  molecules 
replace  KI  molecules  in  the  nanocrystal  cation  of  the  type 
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with  n  =  0,  1,  ....  14.  We  found  that  the  observed  intensity  distribution  is 
largely  determined  by  the  number  of  configurations  (isomers)  that 
different  mixed  nanocrystals  can  have.  A  simple  kinetic  model  shows  that 
this  situation  can  arise  if  the  stability  (internal  energy)  of  the  different 
mixed  nanocrystals  is  independent  of  their  exact  composition,  i.e. 
independent  of  the  value  of  n  within  these  mixed  nanocrystals.  A 
theoretical  calculation  of  the  internal  energy  of  these  nanocrystals  using 
Monte  Carlo  technique  indeed  support  this  assumption. 


II.  Experimental 

The  mass  spectra  of  the  mixed  alkali  halide  nanocrystals  of  the  type 
(RbjjK (n=0-14)  were  obtained  by  the  use  of  a  VG  Analytical 

ZAB-SE*  mass  spectrometer  (reverse  geometry)  fitted  with  a  fast  atom  (Xe) 
bombardment  gun.  The  details  of  the  experimental  have  been  reported 
previously^ ^  Briefly,  a  non-saturated  solution  of  equal  molar  ratio  of  each 
of  Rbl  and  KI  was  prepared  in  distilled  water  and  was  placed  on  the  probe 
tip.  The  water  was  dried  off  leaving  mixed-crystals  on  the  tip^^.  Cluster- 
ions  were  prepared  by  the  bombardment  of  the  mixed-crystals  by  an  8 
keV  Xe  atom  gun.  An  acceleration  voltage  of  8  keV  was  applied  to  the 
cations  to  extract  them  from  the  source.  The  accelerating  ions  were 
analyzed  with  the  magnetic  sector  in  conjunction  with  the  electric  sector, 
and  detected  by  a  conversion  dynode  followed  by  an  electron  multiplier^ 

The  relative  intensity  of  a  certain  mixed  nanocrystal  in  the 
distribution  is  determined  by  the  ratio  of  its  intensity  to  the  sum  of  the 
intensities  of  the  fifteen  peaks  observed  for  all  of  the  nanocrystals  from 
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n  =  0  to  n  =  14.  The  pressure  in  the  acceleration  region  was  kept  at 

~  5  XIO'”  Torr  and  in  the  analyzer  at  -5  XIO’  Torr.  We  collected  three 

sets  of  200  scans  on  three  different  days  to  ensure  reproducibility.  A 

typical  spectrum  is  given  in  Figure  1,  where  the  peaks  corresponding 
to  n  =  0  to  n=14  of  the  nanocrystals  are  marked. 


III.  Results 

We  used  fast  atom  bombardment  mass  spectrometry  to  generate  the 
mixed  cluster-ions  of  the  type  [^^n^l4-n^l3]^  (n=0-14).  A  mass  spectrum 

of  these  cluster-ions  is  shown  in  figure  1.  Using  this  spectrum,  we 
calculated  the  observed  intensity  distribution  for  these  clusters  by 
calculating  fraction  (Fk)  of  ions  observed  for  the  kth  mixed  cluster-ion, 
having  an  intensity  Ik  in  the  mass  spectrum  by  use  of  the  following 
equation: 

Fk  =  (Ik/  X  In  >  X  100  .  (1) 

n=0 

A  'bimodal'  distribution  for  the  probability  distribution  of  these 
cluster-ions  centered  near  n  =  7  was  obtained  and  is  shown  in  Fig.  2. 
These  results  suggests  that  more  nanocrystal  ions  are  detected  with  near 
comparable  number  of  Rbl  and  KI  molecules.  These  are  also  the 
nanocrystals  with  the  maximum  number  of  configurations  (isomers)  that 

one  can  write  for  the  nanocrystals  made  of  3X3X3  lattices  of  atomic 
species.  For  a  nanocrystal  [RbnI^14-nIl3]^’  number  Wj^.  is  given  by  : 

W^=  (14!  /  k!  (14-k)!  )  .  (2) 
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The  fraction  Pj^  is  the  ratio  of  the  number  of  possible  configurations 
(isomers)  for  any  nanocrystal  k  with  k  Rb  to  the  total  number  of  possible 
configurations  for  all  of  the  nanocrystals,  i.e. 

Pk  =(Wk/  XW„)  .  (3) 

n=0 

Figure  3  shows  a  plot  of  the  fraction  Pj^  as  a  function  of  n.  The  maximum  of 
this  distribution  is  at  n  =  7. 

Below,  we  discuss  a  kinetic  model  that  might  be  responsible  for  the 
formation  of  these  mixed  nanocrystals  by  the  sputtering  technique.  We 
then  find  out  the  conditions  under  which  the  kinetic  model  lead  to  a 
distribution  that  is  controlled  by  the  number  of  isomers  that  each  of  these 
mixed  nanocrystals  can  have. 


IV.  Discussion 


A.  Simple  Kinetic  Model  for  Mixed  Nanocrystal  Formation 

Let  us  assume  that  [Ri>n^l4-n^l3]^  nanocrystals  are  formed  rapidly 

from  the  collisions  between  the  atoms  and  ions  which  are  produced 
initially  through  sputtering  of  the  solid  surface.  The  rapid  formation  of 
these  nanocrystals  is  followed  by  cooling  processes  resulting  from 
the  different  channels  of  evaporation  (  unimolecular  decompositions  ).  In 
this  simplified  model,  we  may  write  the  number  of  ions  (  Nk  )  observed 
for  a  certain  cluster  k  as 
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--  formed  ..  decomposed 

Nk  =  -  Nfc 


(4) 


A  certain  cluster  k  can  be  formed  by  nucleation  resulting  from  the 
many  sticky  collisions  between  a  certain  cation  and  the  dipolar  KI  and  Rbl 
molecules.  It  can  also  be  formed  from  the  evaporative  cooling  of  a  less 
stable  nanocrystal  that  has  one  or  two  molecules  of  KI  and/or  Rbl  beyond 
the  magic  number  (  e.  g.  where  n  =  0- 

14). 


B.  Formation  of  Nanocrystals  Through  Nucleation 

The  experiments  of  Ens^^  et.  al.  have  clearly  shown  that  the  intensity 
of  mass  peaks  decrease  exponentially  with  size  if  the  clusters  are 
detected  immediately  after  their  formation  without  giving  them  time  to 
undergo  evaporation.  This  suggests  that  sputtering  atomizes  the  solid, 
and  then  the  atoms  nucleate  to  form  the  different  nanocrystals.  The 
larger  the  number  of  collisions  required  for  the  formation  of  a  nanocrystal, 
the  smaller  is  its  probability  of  formation.  When  nanocrystals  undergo 
evaporative  cooling  prior  to  their  detection,  the  naoncrystals  with  one  or 
two  alkali  halide  molecules  more  than  the  ones  with  magic  numbers 
evaporate^^’  their  excess  molecules  to  attain  the  stable  structure  of  the 
magic  ones. 


It  is  reasonable  to  assume  that  the  nanocrystals  under  study 
are  formed  by  the  heterogeneous  nucleation  in  which  a  metal  cation(  K'*'  or 
Rb'*')  seeds  the  nucleation  of  nanocrystal  by  attracting  KI  and  Rbl 
molecules.  This  would  require  a  certain  number  of  collisions  ,  ,  to  form  a 
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nanocrystal  of  a  fixed  composition  and  arrangement  (  i.e.  a  certain 

isomer).  If  one  assumes  that  a  nanocrystal  of  a  certain  composition(  a  Hxed 
value  of  n,  say  k)  requires  number  of  collisions,  then  the  number  of 


nanocrystal  ions  formed  from  the  nucleation  process,  is  written  as 

j4^form  ^  N^isom  -  (5) 


If  the  collisions  are  all  sticky,  one  can  assume  that  the  number  of 

collisions  required  to  form  any  one  of  the  configurations  of  any  of  the 
fifteen  nanocrystals  is  the  same  (  i.e.  is  the  same  for  any  isomer 

having  a  value  of  n).  This  allows  us  to  write  equation  (5)  as 


Nj^form^  K 


where  the  constant  K  is  equal  to  KV^k^k. 

If  the  detection  time  is  shorter  than  evaporation  time  from  these 
nanocrystals  as  in  the  Ens^^  et.  al.  experiment,  and  if  the  sensitivity  of  the 
detector  is  similar  for  all  the  mixed  nanocrystal,  then  the  predicted 
relative  peak  intensity  is  given  by 


(N^‘ 


14 

I  ) 

n=0 


As  shown  in  eq.  (3),  the  number  of  configurations  (isomers)  that  can 
be  written  for  the  various  compositions  should  have  a  maximum  value 
for  n=7.  This  is  in  reasonable  agreement  with  the  observation  in  Fig.  2. 
However,  under  our  experimental  condition,  it  was  shown  previously  that 
evaporation  indeed  takes  place.  Thus  we  need  to  include  the  effect  of  such 
processes  on  the  nanocrystal  distribution.  This  is  done  in  section  C. 
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C.  Inclusion  of  Evaporative  Cooling 


Incorporation  of  evaporative  cooling  changes  the  kinetic  scheme  by 
introducing  additional  sources  of  formation  as  well  as  a  decay  route  for 
each  nanocrystal.  These  processes  give  rise  to  the  observed  large  relative 

intensities  for  nanocrystals  with  magic  numbers.  Our  previous  studies^  on 
[CSjjIjj  j]'*'  ,  carried  out  on  the  same  instrument  we  are  presently  using, 
have  shown  that  larger  clusters  undergo  evaporation  of  Csl  and  (Csl)2  to 

form  nanocrystals  with  magic  numbers.  'Hot'  nanocrystals  with  magic 
numbers  undergo  losses  of  monomers  and  dimers  but  with  much  reduced 
probability. 

We  assume  a  dominance  of  monomer  alkali  halide  evaporation  as 
well  as  the  steady  state  approximation  in  the  number  of  any  one  of  the 
nanocrystals,  e.  g.  [Rb|jKj4.|jIj3]''' which  we  call  Nj^  13-  This  nanocrystal  is 
formed  by  nucleation  with  a  rate  of  and  from  the  evaporation  of 

a  KI  molecule  from  the  nanocrystal  with  a  rate  of  Pi4Nj^  ,  where 

p  j  4  is  the  rate  constant  and  j  ^  is  the  steady  state  number  of  the 

nanocrystal.  The  [Rbj^Kj4.j^Ij3]'''  nanocrystal  would  give  the 
[Rbj^K|3_jjIj2]'*’  nanocrystal  upon  the  evaporation  of  a  KI  molecule. 


With  the  above  kinetic  scheme,  one  obtains  for  the  steady  state 
number  of  the  M^Ij3''' nanocrystal,  Nj^  ^3  ,  the  following  expression; 


14 

IK^  N 
j=13 


k.j 


isom 


Pk,13 


(8) 


Due  to  the  fact  that  Kj  falls  rapidly  in  value  with  j,  as  seen  from  eq.  (5), 
j=13  will  be  the  dominant  term  in  the  sum  of  eq.  (8). 
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If  the  sensitivity  of  the  detector  is  the  same  for  all  the  mixed 
nanocrystals,  the  predicted  intensity  of  the  peak  corresponding  to  the  k,13 
nanocrystal  is  given  by 

^k.l3  “  13  ^  Pk.l3  . 

The  predicted  intensity  from  equation  (9)  would  account  for  the 
observation  in  Fig.  2,  if  P|j  J3  is  independent  of  k,  implying  that  the  rate 

constant  of  evaporation  from  any  of  the  magic  mixed  nanocrystals  is 
insensitive  to  their  composition. 

For  a  certain  channel  of  evaporation  ,  e.g.  the  loss  of  KI,  the  activation 
energy,  which  determines  the  value  of  the  rate  constant  (p),  is  determined 
by  the  difference  between  the  internal  energies  of  the  and 

the  sum  of  the  energies  of  KI  +  j^Ij2]^-  This  is  only  true  in  the 

absence  of  a  reverse  activation  energy,  which  is  found  to  be  the  case  for 
the  evaporation  of  Csl  from  [Cs^Ij^j. j]'*'  clusters^.  For  Pj^  13  to  be 

independent  of  k,  one  is  thus  forced  to  conclude  that  the  internal 
energies  of  the  nanocrystals  are  similar  and  insensitive  to 

their  composition,  i.e.  to  the  value  of  n.  In  order  to  test  this  conclusion,  we 
have  carried  out  in  the  next  section  a  Monte  Carlo  simulation  of  the 
internal  energies  of  these  clusters  for  n=  0  to  14. 


V.  A.  Computational  Method 

The  Metropolis  Monte  Carlo  method"*^  was  used  to  minimize  the  total 
energy  of  mixed  nanocrystals  through  optimization  of  their  geometry.  We 
used  the  polarizable-ion  model^^’  which  includes  the  deformation 
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of  the  electron  density  as  atoms  come  close  to  one  another.  Even  though 
there  is  a  substantial  free  surface  present  in  these  clusters,  the  covalent 
bonding  forces  are  absent,  and  these  alkali  halides  are  ionically  bound^^. 
Thus,  we  can  use  very  simple  model  potentials  to  study  their  binding 
energies  theoretically^^. 

The  binding  energies  of  the  mixed  nanocrystals  were  calculated  by 
minimizing  the  total  potential  energy  written  as  sum  of  the  pairwise 
coulomb  interaction,  the  short-range  Born-Mayer  repulsive  interaction  due 
to  the  overlap  of  electron  clouds  on  two  nearest  neighbors,  charge-dipole 
interactions,  dipole-dipole  interactions,  and  self-energy  of  the  dipoles^^.  In 
the  Born-Mayer  repulsive  term,  an  effective  radius  is  used  which  accounts 
for  the  deformation  of  the  electron  shells  at  small  distances^  The 
parameters  used  in  our  calculations,  e.g.  polarizabilities,  and  effective  radii, 
were  taken  from  Welch^'*  et.  al.  and  Born^*.  In  the  Born-Mayer  potential, 
the  interactions  for  the  like-charges  have  not  been  considered.  These 
interactions  are  negligible  since  the  overlap  of  the  repulsive  electron 
densities  among  them  is  minimized  due  to  their  relatively  large 
separations. 

In  the  calculation,  we  start  with  the  bond  distances  of  (KI)2 

and/or  (Rbl)2  as  the  initial  bond  distances  between  the  ions  in  our 
nanocrystals^^’  The  nanocrystal  configuration  is  allowed  to  relax  by 
randomly  moving  the  three  independent  coordinates  of  each  atom  in  the 
nanocrystal.  After  each  move,  a  change  in  the  total  energy  (  AE)  between 
the  initial  and  final  structure  is  calculated.  If  AE  is  negative  or  zero,  the 
new  structure  is  accepted,  otherwise,  it  is  accepted  with  a  probability  of 
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exp(-AE/kT)^^’^^.  The  temperature  is  lowered  gradually  for  the  calculation 
of  the  ground  state  energy.  About  200000  (  or  more)  sweeps  are 
performed  in  order  to  increase  the  probability  of  reaching  a  global 
minimum^'*.  The  computer  code  for  our  calculation  was  modified  from  that 
of  Li^^.  The  results  of  our  calculations  for  the  binding  energies  of  all  the 
different  nanocrystals  are  given  in  Table  1. 


V.  B.  Computational  Results 

The  results  of  our  calculation  for  the  total  binding  energies  are  given 
in  Table  1.  A  slight  monotonic  increase  in  the  calculated  values  of  the 
binding  energies  is  obtained  as  the  number  of  Rb  ions  in  the  mixed 
nanocrystal  increases.  However,  the  difference  in  the  energies  is  smaller 
than  the  inherent  error  in  our  calculations.  Thus,  within  the  limit  of  the 
errors  in  our  calculations,  we  can  conclude  that  the  internal  energies  of 
the  mixed  nanocrystals  of  the  type  [RbAKH-nllS]"*"  (n=0-14)  are  rather 
insensitive  to  their  composition. 

The  vibrational  frequencies  of  the  monomeric  KI  is  larger  than  that 
for  Rbl  (186.53  cm'l,  and  138.511  cm-1,  respectively)'*'*  in  the  gas  phase. 
The  effect  of  this  would  come  in  the  zero  point  energy  contributions  to  the 
binding  energies  of  the  nanocrystals  of  [RbnKl4-nIl3]‘'’  (n=0-14).  This  is 
expected  to  decrease  the  calculated  energy  difference  between 
nanocrystals. 

Furthermore,  in  our  calculation  of  the  binding  energies,  we 
considered  only  a  small  number  of  the  possible  isomers  for  each  of  the 
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[RbnKl4-nIl3]‘’’  (n=0-14)  nanocrystals.  There  are  many  isomers  for  any 

one  of  the  nanocrystals  and  many  of  these  isomers  would  be  present 

significantly  at  the  high  temperatures  of  our  experiment,  which  is  -  750  K. 

All  these  factors  suggest  that  the  small  calculated  difference  in  the  internal 

energies  of  the  nanocrystals  is  not  significant.  One  can  then  conclude  that 

these  mixed  nanocrystals  have  comparable  energies.  This  conclusion  is  not 

22 

at  all  surprising  in  the  light  of  the  fact  that  the  lattice  energy  of  bulk  KI 
(-I49.9kcaI/mol)  is  comparable  to  that  of  Rbl(  -144.9  kcal/mol). 
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Table  1  The  calculated  binding  energies  of  the  [RbjjKj4_j^Ij3]+  (n=0-14) 
nanocrystals  as  obtained  through  Monte  Carlo  Simulation  Annealing. 


n  in  the  [RbnKi4.„Ii3]+ 

0 

1 

2 

3 

4 
6 

7 

8 

9 

10 
1 1 
12 

13 

14 

Average  Energy 
*  average  of  more  than  one 


Binding  Energy  (eV) 

-80.52±1.8 
-80.65+3.62* 
-80.54+3.60* 
-80.10±3.58* 
-79.76±3.57 
-79.78±3.55* 
-79.53±3.55* 
-78.33±3.51* 
-78.34±3.51 
-78.26±3.50 
-78.06±3.44 
-78.37±3.45* 
-78.59±3.60* 
-77.96±3.50* 
-79,23±3.50 
isomer. 
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Figure  Captions 


Figure  1  The  observed  mass  spectrum  resulting  from  the  sputtering  of  a 
solid  mixture  of  KI  and  RbL  Marked  are  the  peaks  corresponding  to  the 
mixed  magic  nanocrystals  [RbnKi4-nIl3]‘*'  for  n  =  0,  1, ....  ,  14  . 


Figure  2  The  calculated  relative  intensity  of  kth  mixed  nanocrystal 
(relative  to  all  the  nanocrystals  formed)  as  calculated  from  the  observed 
intensities  of  the  fifteen  peaks  of  the  nanocrystals  of  [RbnKi4-nIl3]‘^  for 
n=  0-14- 


Figure  3  The  predicted  dependence  of  the  relative  number  of  isomers 
of  the  mixed  nanocrystals  on  their  composition(i.e.  the  number  of  KI  and 
Rbl  molecules)  .  The  mixed  nanocrystal  having  equal  number  of  KI  and  Rbl 
(ns7)  has  the  largest  relative  number  of  isomers. 
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